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Different characteristic of matter influencing the evolution of the Universe has been sim- 
ulated by means of a nonlinear spinor field. Exploiting the spinor description of perfect fluid 
and dark energy evolution of the Universe given by an anisotropic Bianchi type- VI, VIo, V, 
HI, I or isotropic Friedmann-Robertson-Walker (FRW) one has been studied. It is shown that 
\ due to some restrictions on metric functions, initial anisotropy in the models Bianchi type- 

5_i \ VI, VIo, V and III does not die away, while the anisotropic Bianchi type-I models evolves 

into the isotropic one. 
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I. INTRODUCTION 

Being related to almost all stable elementary particles such as proton, electron and neutrino, 
spinor field, especially Dirac spin-1 /2 play a principal role at the microlevel. However, in cos- 
mology, the role of spinor field was generally considered to be restricted. Only recently, after 
some remarkable works by different authors (e.g., Henneaux, 1980; Saha, 2001a, 2004a, 2004b, 
2006c, 2006d; Armendariz-Picon and Greene, 2003; Ribas et al, 2005; Souza and Kremer, 2008), 
showing the important role that spinor fields play on the evolution of the Universe, the situation 
began to change. This change of attitude is directly related to some fundamental questions of 
modem cosmology: (i) problem of initial singularity; (ii) problem of isotropization and (iii) late 
time acceleration of the Universe. 

(i) Problem of initial singularity: One of the problems of modern cosmology is the presence 
of initial singularity, which means the finiteness of time. The main purpose of introducing a 
nonlinear term in the spinor field Lagrangian is to study the possibility of the elimination of initial 
singularity. In a number of papers it was shown that the introduction of nonlinear spinor field 
into the system indeed gives rise to singularity-free models of the Universe (Saha, 2001a, 2001b, 
2004a), (Saha and Shikin, 1997a, 1997b). 

(ii) problem of isotropization: Although the Universe seems homogenous and isotropic at 
present, it does not necessarily mean that it is also suitable for a description of the early stages 
of the development of the Universe and there are no observational data guaranteeing the isotropy 
in the era prior to the recombination. In fact, there are theoretical arguments that support the 
existence of an anisotropic phase that approaches an isotropic one (Misner 1968). The observations 
from Cosmic Background Explorer's differential radiometer have detected and measured cosmic 
microwave background anisotropics in different angular scales. These anisotropics are supposed to 
hide in their fold the entire history of cosmic evolution dating back to the recombination era and are 
being considered as indicative of the geometry and the content of the universe. More about cosmic 
microwave background anisotropy is expected to be uncovered by the investigations of microwave 
anisotropy probe. There is widespread consensus among the cosmologists that cosmic microwave 
background anisotropics in small angular scales have the key to the formation of discrete structure. 
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It was found that the introduction of nonlinear spinor field accelerates the isotropization process 
of the initially anisotropic Universe (Saha, 2001a, 2004a, 2006c). 

(iii) late time acceleration of the Universe: Detection and further experimental reconfirmation 
of current cosmic acceleration pose to cosmology a fundamental task of identifying and revealing 
the cause of such phenomenon. This fact can be reconciled with the theory if one assumes that the 
Universe id mostly filled with so-called dark energy. This form of matter (energy) is not observable 
in laboratory and it does not interact with electromagnetic radiation. These facts played decisive 
role in naming this object. In contrast to dark matter, dark energy is uniformly distributed over the 
space, does not intertwine under the influence of gravity in all scales and it has a strong negative 
pressure of the order of energy density. Based on these properties, cosmologists have suggested a 
number of dark energy models, those are able to explain the current accelerated phase of expansion 
of the Universe. In this connection a series of papers appeared recently in the literature, where a 
spinor field was considered as an alternative model for dark energy (Ribas et al, 2005; Saha, 
2006d, 2006e, 2007). 

It should be noted that most of the works mentioned above were carried out within the scope 
of Bianchi type-I cosmological model. Results obtained using a spinor field as a source of 
Bianchi type-I cosmological field can be summed up as follows: A suitable choice of spinor field 
nonlinearity 

(i) accelerates the isotropization process (Saha, 2001a, 2004a, 2006c); 

(ii) gives rise to a singularity-free Universe (Saha, 2001a, 2004a, 2004b, 2006c); 

(iii) generates late time acceleration (Ribas etal, 2005; Saha, 2006d; Souza and Kremer, 2008). 
Given the role that spinor field can play in the evolution of the Universe, question that naturally 

pops up is, if the spinor field can redraw the picture of evolution caused by perfect fluid and 
dark energy, is it possible to simulate perfect fluid and dark energy by means of a spinor field? 
Affirmative answer to this question was given in the a number of papers (Krechet et al, 2008; 
Saha, 2010a, 2010b). In those papers the authors have shown that different types of perfect fluid 
and dark energy can be described by nonlinear spinor field. In (Saha, 2010a) we used two types 
of nonlinearity, one occurs as a result of self-action and the other resulted from the interaction 
between the spinor and scalar field. It was shown that the case with induced nonlinearity is the 
partial one and can be derived from the case with self-action. In (Saha, 2010b, 2011) we give the 
description of generalized Chaplygin gas and modified quintessence in terms of spinor field and 
study the evolution of the Universe filled with nonlinear spinor field within the scope of a Bianchi 
type-I and FRW cosmological model. The purpose of this paper is to extend that study within the 
framework of other Bianchi models. 



II. SIMULATION OF PERFECT FLUID WITH NONLINEAR SPINOR FIELD 

Nonlinear quantum Dirac fields were used by Heisenberg (1953, 1957) in his ambitious unified 
theory of elementary particles. They are presently the object of renewed interest since the widely 
known paper by Gross and Neveu (1974). A nonlinear spinor field, suggested by the symmetric 
coupling between nucleons, muons, and leptons, has been investigated by Finkelstein et al. (1951) 
in the classical approximation. 

In this section we simulate different types of perfect fluid and dark energy by means of a 
nonlinear spinor field. 



A. Spinor field Lagrangian 

For a spinor field y, the symmetry between y and \[r appears to demand that one should choose 
the symmetrized Lagrangian (Kibble, 1961). Keeping this in mind we choose the spinor field 
Lagrangian as (Saha, 2001a): 



(2.1) 
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where the nonlinear term F describes the self-action of a spinor field and can be presented as 
some arbitrary functions of invariant generated from the real bilinear forms of a spinor field. For 
simplicity we consider the case when F = F(S) with S = ty/ty/. Here is the covariant derivative 
of spinor field: 

v MV = ^r-r M v, v M v/=^ + ^r M , (2.2) 

with being the spinor affine connection. In (12.11 ) y's are the Dirac matrices in curve space-time 
and obey the following algebra 

y n r v + f r n = 2 g nv (2 . 3) 
and are connected with the flat space-time Dirac matrices f in the following way 

8pv(x) = eJ(x)4(*)T?a&, 7m W = e % (x)?a, (2-4) 

where r\ a b = diag(l, — 1, — 1, — 1) and is a set of tetrad 4- vectors. The spinor affine connec- 
tion matrices (x) are uniquely determined up to an additive multiple of the unit matrix by the 
equation 

v M 7v = -r^r P -r M 7v + r v r M = o, (2.5) 

with the solution 

T^ = \yafd^ ] - l -y p fY%, (2.6) 
Varying (12.11 ) with respect to y/tyA) one finds the spinor field equations: 

dF 

ifV^-m sp f+-y/ = 0, (2.7a) 
ii dF 

N^f+m sp xfr-—\fr = 0, (2.7b) 
The energy-momentum tensor of the spinor field is given by 

T£ = ^g pv (n V v W + W7v V M W - V M \fry v y - V v yty, yrj - 5/?L sp (2.8) 
where L sp in account of spinor field equations (|2.7al ) and (I2.7bl ) takes the form 

L SV = -S^ + F(S). (2.9) 

We consider the case when the spinor field depends on t only. In this case for the components 
of energy-momentum tensor we find 

r ° = m sp 5-F, (2.10a) 
T{ = T? = T 3 3 = S^-F. (2.10b) 

A detailed study of nonlinear spinor field was carried out in Saha (2001a, 2004a, 2006c). In what 
follows, exploiting the equation of states we find the concrete form of F which describes various 
types of perfect fluid and dark energy. 
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B. perfect fluid with a barotropic equation of state 

First of all let us note that one of the simplest and popular model of the Universe is a homo- 
geneous and isotropic one filled with a perfect fluid with the energy density e = r o ° and pressure 
p = — Tj 1 = — T 2 2 = — T 3 3 obeying the barotropic equation of state 

p = We, (2.11) 

where W is a constant. Depending on the value of W (|2.1 II) describes perfect fluid from phantom 
to ekpyrotic matter, namely 

W = 0, (dust), (2.12a) 

W = 1/3, (radiation), (2.12b) 

W E (1/3,1), (hard Universe), (2.12c) 

W = 1, (stiff matter), (2.12d) 

W G ( — 1/3,-1), (quintessence), (2.12e) 

W = — 1, (cosmological constant), (2.12f) 

W < —1, (phantom matter) , (2.12g) 

W > 1, (ekpyrotic matter). (2.12h) 

The barotropic mo del o f perfect fluid is used to study the evolution of the Universe. Most re- 
cently the relation (|2.1 II) is exploited to generate a quintessence in order to explain the accelerated 
expansion of the Universe (Saha, 2005; Zlat ev 199 9). 

In order to describe the matter given by (|2.12t with a spinor field let us now substitute £ and 
p wi th r n ° and — 7} , respectively. Thus, inserting e = Tq and p = —T^ from (|2. 10a!) and (|2.10bl) 
into (12.111) we find 

S^-(\+W)F + m sp WS = 0, (2.13) 
with the solution (Saha, 2010a, 2010b, 2011) 

F = XS l+w +m sp S. (2.14) 

Here A is an integration constant. Taking into account that the energy density should be non- 
negative, we conclude that A is a negative constant, we write the components of the energy mo- 
mentum tensor 

r ° = vS l+w , (2.15a) 

r 1 1 = r 2 2 = r 3 3 = -vws l+w , (2.15b) 

where X = — V, with v being a positive constant. As one sees, the energy density e = r o ° is always 

positive, while the pressure p = —T* = vWS l+w is positive for W > 0, i.e., for usual fluid and 
negative for W < 0, i.e. for dark energy. 

In account of it the spinor field Lagrangian now reads 



L -L 



vS l+w , (2.16) 



Thus a massless spinor field with the Lagrangian (12.161 ) describes perfect fluid from phantom to 
ekpyrotic matter. Here the constant of integration v can be viewed as constant of self-coupling. A 
detailed analysis of this study was given in Krechet (2008). 
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C. Chaplygin gas 

An alternative model for the dark energy density was used by Kamenshchik et al. (2001), 
where the authors suggested the use of some perfect fluid but obeying "exotic" equation of state. 
This type of matter is known as Chaplygin gas. The fate of density perturbations in a Universe 
dominated by the Chaplygin gas, which exhibit negative pressure was studied by Fabris et al. 
(2002). Model with Chaplygin gas was also studied in the Refs. (Dev et al, 2003; Bento et al, 
2002). 

Let us now generate a Chaplygin gas by means of a spinor field. A Chaplygin gas is usually 
described by a equation of state 

p = -A/e r . (2.17) 
Then in case of a massless spinor field for F one finds 

(-F)rd(-F) _ dS 

(_ F )l+y_ A S ' 1 ' 

with the solution (Saha, 2010a, 2010b, 201 1) 

-F=(A + A5 1+ ^) 1/(1+r) . (2.19) 
On account of this for the components of energy momentum tensor we find 

r ° = (A + A5 1+ ^) 1/(1+r) , (2.20a) 
T} = Ti = T 3 3 = A/(A + XS l+ V) Y/{l+Y) . (2.20b) 

As was expected, we again get positive energy density and negative pressure. Analogical results 
were obtained in (Cai et al, 2008 ). 

Thus the spinor field Lagrangian corresponding to a Chaplygin gas reads 



L -L 



(A + XS l+ vy /[l+Y) . (2.21) 



Setting 7 = 1 we find the result obtained in Saha (2010a). 



D. Modified quintessence 

Finally, we simulate modified quintessence with a nonlinear spinor field. It should be noted that 
one of the problems that face models with dark energy is that of" eternal acceleration. One of the 
possible way to avoid this problem is to introduce a negative A-term together with a quintessence, 
which gives rise to a oscillatory mode of expansion (Cardenas et. al, 2003, Saha, 2006a) In order 
to get rid of that problem quintessence with a modified equation of state was proposed which is 
given by (Saha, 2006b) 

p = W(e-e a ), We (-1,0), (2.22) 

Here e cr some critical energy density. Setting e cr = one obtains ordinary quintessence. It is well 
known that as the Universe expands the (dark) energy density decreases. As a result, being a linear 
negative function of energy density, the corresponding pressure begins to increase. In case of an 
ordinary quintessence the pressure is always negative, but for a modified quintessence as soon as 
£ becomes less than the critical one, the pressure becomes positive. 
Inserting e = r o ° and p = —T\ into (12.221 ) we find 

F = -vS l+w +m sp 5- — -£ cr , (2.23) 
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with rj being a positive constant. On account of this for the components of energy momentum 
tensor we find 



vS 



i+w 



+ 



W 



l+W 



-CI: 



-vWS 



i+w 



+ 



W 



\+w 



(2.24a) 
(2.24b) 



Lagrangian for spinor field describing perfect fluid and modified quintessence can be written 
in the following way (Saha, 201 1) 



vS 



i+w 



W 



l+W 



(2.25) 



One can easily verify, in case of £ cr = (12.251) corresponds to a perfect fluid, while nontrivial e cr 
with WE (—1,0) generates modified quintessence. It should be noted that the restriction W > — 1 
is very important, as it does not allow the system to cross over phantom divide barrier. 

We see that a nonlinear spinor field with specific type of nonlinearity can substitute perfect 
fluid and dark energy, thus give rise to a variety of evolution scenario of the Universe. 



III. COSMOLOGICAL MODELS WITH A SPINOR FIELD 

In the previous section we showed that the perfect fluid and the dark energy can be simulated by 
a nonlinear spinor field. In the section II the nonlinearity was the subject to self-action. In (Saha, 
2010a) we have also considered the case when the nonlinearity was induced by a scalar field. It 
was also shown the in our context the results for induced nonlinearity is some special cases those 
of self-interaction. Taking it into mi nd w e study the evolution an Universe filled with a nonli near 
spinor field given by the Lagrangian (|2.1I ). with the nonlinear term F is given by (|2.14l) . (|2.19l) and 



A. Bianchi type anisotropic cosmological model 

Let us study the evolution of an anisotropic Bianchi type cosmological model filled with spinor 
field. In this report we consider Bianchi type- VI, VIo, V, HI, I and FRW models. 
We choose the Bianchi type- VI cosmological model in the form (Saha, 2004b) 

ds 2 = dt 2 - of e - 2mz dx 2 - a\e 2nz dy 2 - ajdz 2 , (3.1) 

with a\, ci2, as being the functi on of t only. A suitable choice of m, n as well as the metric functions 
a\, ci2, CZ3 in the BVI given by (13.11) evokes the following Bianchi-type universes: 

• for m = n the BVI metric transforms to a Bianchi-type VIo (BVIo) one, i.e., m = n, BVI =^> 
BVIo G open FRW with the line elements 

ds 2 = dt 2 - a\e- 2mz dx 2 - a\e 2mz dy 2 - a\ dz 2 ; (3.2) 

• for m = —n the BVI metric transforms to a Bianchi-type V (BV) one, i.e., m = n, BVI =^ 
BV G open FRW with the line elements 



ds 2 = dt 2 - a\e 2mz dx 2 - a\e 2mz dy 2 - a\ dz 2 ; 



(3.3) 
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• for n = the BVI metric transforms to a Bianchi-type III (Bill) one, i.e., n = 0, BVI =>- 
Bill with the line elements 

ds 2 = dt 2 -a\e- 2mz dx 2 -a 2 2 dy 2 -aldz 2 \ (3.4) 

• for m = n = the BVI metric transforms to a Bianchi-type I (BI) one, i.e., m = n = 0, BVI 
=>• BI with the line elements 

ds 2 = dt 2 — a\dx 2 — a 2 dy 2 — a 2 dz 2 ; (3.5) 

• for m = n = and equal scale factor in all three directions the BVI metric transforms to 
a Friedmann-Robertson-Walker (FRW) universe, i.e., m = n = and a = b = c, BVI =>- 
FRW with the line elements. 

ds 2 = dt 2 -a 2 (dx 2 + dy 2 + dz 2 ). (3.6) 

Let us go back to the metric Bianchi-type VI. The metric (13.11 ) possesses following nontrivial 
co variant and contravariant components: 

goo = 1, gn = a\e- 2mz , g 22 = aje 2nz , g 33 = a% 

a°0 =1 „H = \ g 22 = 1 a33 = _L 

8 ' 8 a\e-^ 8 a\e 2 ^ 8 a\ 

The nontrivial Christoffel symbols for (13.11 ) are 

x _ fl'i 2 _ d 2 r 3 _ a '3 
a\ 0-2 d 3 

T° n = a x d x e- 2m \ I% 2 = a 2 d 2 e 2nz , T%=a 3 d 3l (3.7) 

-.2 m „2 



T 1 - m r 2 -n T 3 _ ma l --2mz r 3 _ _ na 2 2nz 
1 31 — m ' 1 32 — n 7 1 11 — — 2~ e ' 1 22~ T e • 



In view of (12.41) we choose the tetrad as follows: 

4 0) = 1, e[ l) =a ie - mz , ef=a 2 e nz , ef = a 3 . (3.8) 

From 

J» = ef%. (3.9) 

one now finds 

7o = 7o, 7i =aie~ mz % y 2 = a 2 e nz % y 3 = a 3 f 3 . (3.10) 
Taking into account that in our case 

f = % f = -fi, f = ~% ? = -% 

one also finds 



e mz e~ nz 1 
7° = 7°, f = — y\ y 2 = f, y i = —f. 

<x\ a 2 a 3 



(3.11) 
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Now we are ready to compute spinor affine connections using (12.61) which in our particular case 
gives 



r = 


1 v ^ (a) 


-^7p7 v rg v =o, 


(3.12a) 


n = 


-\r P fr p lv 




(3.12b) 


r 2 = 


-\r P fr p 2v 


2 V a 3 ) 


(3.12c) 


r 3 = 




-\r P fr p 3v = \a 3 ff. 


(3.1 2d) 



From (13.121 ) one finds 

1/^1 d 2 d 3 \ n. m — n _ % ,„ „ „ 

r M r M = -- -+-+^ V + ^^r 3 - O.13 

2 Vai a2 03/ 2fii3 
Choosing m, n we can thus write the spinor affine connections for other Bianchi type metrics. 



1. Bianchi type-VI anisotropic cosmological model 

Let us now study the evolution of the Universe gi ven b y a Bianchi type-VI cosmological model. 
The Einstein equations corresponding to the metric (|3.1t take the form: 

— + — + — — = kT{, (3.14a) 

a 2 «3 «2 ®3 «3 

«3 . at . d 3 di m 2 2 

1 1 7T = kT 2 , (3.14b) 

a 3 a\ a 3 a\ a 3 

a\ d 2 d\ d 2 , mn 3 
H 1 h— = kT/, (3.14c) 



ai a2 ai a 2 a 3 
di d2 ^2 d 3 d 3 d\ m 2 — mn + n 2 



+ z±z± + z±z±-- — = kT o^ (314d) 
a\ a 2 a 2 a 3 a 3 a\ a 3 

m ^i_ n ^_( m _ n )^3 = kT o, (3 .i4e) 
ci\ a 2 cl 3 

Taking into account that we are dealing with massless spinor field, from (12. lOab and (12.1 Obi ) one 
finds r ° = —F, T} = T 2 = T 3 3 = S% -F and r 3 ° = 0. Then the Eq. (l3T4el) immediately gives 

1 11 M-), Jf\ = const.. (3.15) 

\a 3 / 1 

Let us now define 



a 3 J \a 3 



V = a\a 2 a 3 . (3.16) 
Summation of (|3.14at , (13. 14bb . (I3.14c|) and three times (|3.14dt gives 



V m —mn + n 3k 

V 4 + 



dF ■ 

S IF 

dS 



(3.17) 
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Before solving this equation let us write S in terms of v. The spinor field equations in this case 
read 

nf . V \ m — n_^ .dF „ ,„ „ 



From (13.181 ) one finds 
which gives 



2a 3 

diys) 

dt 



(3.19) 



S = Co = const. (3.20) 

Thus we see that the equation for defining V explicitly depends on a^. Here we assume that 
the expansion is proportional to the eigenvalue a\ of shear tensor a^. In a comoving system of 
reference with U^ 1 = (1, 0, 0, 0) for Bianchi type- VI metric we find 

e = U» = ^ + ^ + ^. (3.21) 

Cl\ U2 «3 

From 

a MV = \{U^+U v , p pP)- l -eP lLV , (3.22) 

we find 

a/ = --(-2- + ^ + ^). (3.23) 
3 v ai a2 ^3 

As one sees, S, and o\ (it is true for o\ and a| as well) do not depend on m and n. Hence, they 
will remain unaltered for other Bianchi models following from (13.11) . Now setting 

ai=tfie, (3.24) 

where q\ is a constant, one finds 

ay = {a 1 a z ) jr \ jV 2 = (1 +3q 1 )/(2-3q 1 ). (3.25) 
Inserting this into (13.151) and using (13.161) one finally finds: 

ai = y^/(i+^ 2 ) 5 (3.26a) 

a 2 = ^- l /( 2n - m )y( m - n --- A 2>n)/(^2+i)(rn-2n) ^ (3.26b) 
a3 — ^l/{2n-m)y{Jim-n)/(,yr 2 +i)(m-2n) ^ (3.26c) 

Let us note that perfect fluid satisfying the barotropic equation of state, as well as Chaplygin 
gas are described by a massless spinor field Lagrangian. Moreover, S, as well as F(S) are the 
functions of V. Therefore, we can rewrite the equations for defining V as 



dF 

S IF 

dS 



V, (3.27) 
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with 



2{n — Jr^m) 1 2qi (m + n) 

qi = {\+jV 2 ){m-2n) + 1 = 3 m-2n ' 



In case of spinor field nonlinearity given by (|2.16l ). Eq. (|3.27l) takes the form 

m 2 -mn + n 2 . 3?cvd +M/ (l - W) w 
V = 2 V qi H - -V 



2/(2w— m) 



(3.28) 



with the solution in quadrature 

dV 



q 3 V42+ l +3KvC^ +w V i - w +Ci 



t+tQ, <?3 



4(m 2 — mn + n z 



(92 + 1K 



7. 1 {In— m) ' 



(3.29) 



Here Ci is some integration constant. 

Let us consider the case when the spinor field is given by the Lagrangian (|2.21l) . The equation 
for V now reads 



m 2_ mn + n 2 3K 

L ^2/(2n-m) + 2 

with the solution 

dV 



(Av l+ y+xc 1 +r ) m+r) + 



AV l+ V 



(AV l +r + ?iC^ r ) 



i+r\r/(i+r) 



(3.30) 



^Ci+q 3 Vi2+ l +3KV(AV l +r + ?LC 1 +r ) 



l+y\ 1/(1+7) 



t+to, Ci = const, to = const. (3.31) 



Inserting y = 1 we come to the result obtained in (Saha, 2005). 

Finally we consider the case with modified quintessence. In this case for V we find 



V = 2 



m 2 — mn + n 2 



Vl2 + — [ 

2/(2n-m) 2 



vd +w '(1 -W)V~ W + 2We CI V /(l +W) 



(3.32) 



with the solution in quadrature 



^ 3 V ?2+1 + 3fc[vC ( 5 +M/ y 1 - H/ +We CI V 2 /(l + W)] +d 



t + tQ. 



(3.33) 



Recalling that in case of modified quintessence W € (— 1, 0). So the model allows cyclic mode of 
expansion, only when q% < 1 . 



2. Bianchi type-VIo anisotropic cosmological model 



Setting m = n from (13.11) we get Bianchi type-VIo cosmological model given by (13.21) . In this 
case equation (I3.14el) gives 

d\ dj 

— - — = 0, (3.34) 
a\ a.2 



with the solution 



a.2 = Jr\a\. 



(3.35) 
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Assuming that a\ °= 0, on account of (13.251 ) in this case we find 

ax = y^/(i+^2) 5 (3.36a) 
a 2 = ^y^/a+^i), (3.36b) 

a 3 = i-y( 1 -^)/(i+^). (3 . 36c) 
JVi 

The equation for V in this case reads 

V = 2m 2 ^ 1 2 V q2 +@(V). (3.37) 

In this case for VIo we find the solution in quadrature analogous to (13.291 ), (13.311 ) and (13.331 ), for 
quintessence, Chaplygin gas and modified quintessence, respectively, with 

2LY 2 -\) 1 A 4m 2 ^ 2 3m 2 ^ 2 



3 *' ^2 + 1 1+3^1 

As we see, in presence of a modified quintessence, BVIo cosmological model allows cyclic 
mode of evolution if q\ < 1/6. In case of q\ > 1/6 we find the model is not bound from above as 
is seen from figure |3] 



3. Bianchi type-V anisotropic cosmological model 

Setting m = — n fr om (13.11 ) we get Bianchi type-V cosmological model given by (13.31) . In this 
case equation (13.14el) gives 

^1 + ^-2^ = 0, (3.38) 

with the solution 

a\a% = ^4ia 2 , (3.39) 
Assuming that a\ 0, on account of (13.251) in this case we find 



a 1 

a 3 



The equation for v in this case reads 



y^ 2 /(i+^2) 5 (340a) 

^l/3 y (2-^)/3(l + ^ 2 ) ? (3 _ 40b) 

-4^/ 3 . (3.40c) 

^ 1/3 



V = 6m 2 ^ 2/ V <?2 + ^(V). (3.41) 

In this case for V we find the solution in quadrature analogous to (13.291) , (13.311) and (13.331 ), for 
quintessence, Chaplygin gas and modified quintessence, respectively, with 

3 <?2 + 1 

We can conclude that in presence of modified quintessence BV model always undergoes a 
cyclic mode of expansion. 
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4. Bianchi type-Ill anisotropic cosmological model 

Setting n = from (13.11 ) we get Bianchi type-VIo cosmological model given by (13.41) . In this 
case equation (I3.14el) gives 

— - — = 0, (3.42) 

with the solution 

a 3 = Jia u (3.43) 
Assuming that o\ °= 0, on account of (13.251 ) in this case we find 

ai = ^ 2 /(i+^) 5 (3 . 44a) 

a 2 = JLy(i-^)/(i+^2) ( 3.44b) 

a 3 = ^y^/(i+-^). (3.44c) 

The equation for V in this case reads 

2 

V = 2—2^ 2 + ®(V). (3.45) 

In this case for V we find the solution in quadrature analogous to (13.291) , (13.311) and (13.331 ), for 
quintessence, Chaplygin gas and modified quintessence, respectively, with 

2jV 2 1 „ 4m 2 6m 2 

<?2 = 7} , : + 1 = -X ~2gl, <?3 



^i + l 3 ^ ( qi +\)jn (2-3^)^2" 

As we see, in presence of a modified quintessence, BVIo cosmological model allows cyclic 
mode of evolution if q\ > —1/3. 

In the figures [TJ|2l|3]and|4]we have plotted the evolution of V for different Bianchi models filled 
with quintessence and modified quintessence. 



5. Bianchi type-I anisotropic cosmological model 

Bianchi type-I (BI) model is the simplest anisotropic cosmological model and gives an excellent 
scope to take into account the initial anisotropy of the Universe. Given the importance of BI 
model to study the effects of initial anisotropy in the evolution of he Universe, we study this 
models in details. Unlike the models considered prev iously, the system of Einstein's equations 
for BI model does not contain off-diagonal component (|3.14el) . explicitly relating metric functions 
between themselves. The system of Einstein equations in this case reads 



a 2 _|_ a 3 a 2 a 3 




(3.46a) 


(12 CI3 &2 a 3 




Q>3 OL\ ®3 Q>\ 


= *T$, 


(3.46b) 


$3 U\ a 3 a\ 






CL\ &2 CL\ CI2 


= kT 3 \ 


(3.46c) 


a\ 0,2 a\ 0:2 






@1 a 2 a 2 a 3 a 3 a \ 


= KT ° 


(3.46d) 


Ol\ Q2 Cl2 a 3 a 3 a l 
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FIG. 1: Evolution of the Universe filled with 
quintessence for different Bianchi models 



FIG. 2: Evolution of the Universe filled with 
modified quintessence for different Bianchi 
models 



Solving the Einstein equation for the metric functions one finds (Saha, 2001a) 

dt^ 3 3 



M = DtV l/3 exp (Xi —), Y\Di ■ = 1 , £ X = 0, 

v J v J i=i (=l 



with Dj and X, being the integration constants. 

The equation for v in this case takes the form (Saha, 2001a) 

v = 9(y). 

In case of (12.161) Eq. (13.481) takes the form 

V = (3/2) kvC 1 +w (1 -W)V~ W 

with the solution in quadrature 

dV 

t + t . 



3kvC 1 + w V 1 -w +Ci 



(3.47) 



(3.48) 



(3.49) 



(3.50) 



Here C\ and to are the integration constants. 

In the Figs. [5] and [6] we have plotted the evolution of the Universe defined by the nonlinear 
spinor field corresponding to perfect fluid and dark energy (Saha, 2010b). 

Let us consider the case when the spinor field is given by the Lagrangian (|2.21l) . The equation 
for V now reads 



V = (3/2)k 



{AV 1+ ? + XC l +7 ) mi+Y) +AV 1+ y{AV l+ V + XC l +Y ) 



i+y\y/(i+y) 



(3.51) 
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2,5 



/ 
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/ 
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/ 

/ BVI 



0,05 



0,15 0,20 t 0,25 




FIG. 3: Evolution of the Universe filled with 
modified quintessence for Bianchi type-VIo 
model 



FIG. 4: Evolution of the Universe filled with 
modified quintessence for Bianchi type-V model 




t 4 




FIG. 5: Evolution of the Universe filled with per- 
fect fluid. 



FIG. 6: Evolution of the Universe filled with 
dark energy. 



with the solution 



dV 



/ c 1 +3 K y(Ayi+r + Ac ( ; +r ) 1/(1+r) 

Inserting y = 1 we come to the result obtained in (Saha, 2005) 



t + to, C\ = const, to = const. 



(3.52) 
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Finally we consider the case with modified quintessence. In this case for V we find 



V = (3/2) K tjCq (1 + W)V W -2We CI V/(l-W) 
with the solution in quadrature 



/ 



dV 



3K[qC} ) - w V l + w -W£ C rV 2 /(l-W)] +Ci 



t + t . 



(3.53) 



(3.54) 



Here C\ and to are the integration constants. Comparing (13.541 ) with those with a negative A-term 
we see that e cr plays the role of a negative cosmological constant. 



Energy density 




Pressure 




FIG. 7: Dynamics of energy density and pres- 
sure for a modified quintessence. 



FIG. 8: Evolution of the Universe filled with a 
modified quintessence. 



In the Fig. [7] we have illustrated the dynamics of energy density and pressure of a modified 
quintessence. In the Fig. [8] the evolution of the Universe defined by the nonlinear spinor field 
corresponding to a modified quintessence has been presented. As one sees, in the case considered, 
acceleration alternates with declaration. In this case the Universe can be either singular (that ends 
in Big Crunch) or regular. 



6. FRW cosmological models with a spinor field 

Since our Universe is almost isotropic at large scale it would be fitting to study evolution of 
the FRW model within the scope of spinor description of matter. The Einstein equations read 
corresponding to the FRW model (13.61 ) reads 



•2 

a a 

2 - + - 

a a L 

„d 2 



kT 



kT ( 



o • 



(3.55a) 
(3.55b) 
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From the spinor field equations in this case we find 



S=— ^, Co = const. 



(3.56) 



In o rder to find the solution that satisfies both (|3.55al) and (|3.55bl) we rewrite (I3.55al ) in view of 
(I3.55bl) in the following form: 

KT/._, _,A (3 _ 57) 



a = -\3Ti-T§\a. 



Further we solve this equation for concrete choice of source field. 

Let u s consi der th e case of perf ect fluid given by the barotropic equation of state. In account of 
(f27T5ab . (l2l5bl) and (13361) . (1537b takes the form " 



kv{\+W)C 1 + w M2+m) 



a = 



(3.58) 



that admits the first integral 



a 2 = ^vC l +w a-^ +3W ^+E u E l = const. 



In Fig. |9] and [10] we plot the evolution of the FRW Universe for different values of W. 



(3.59) 



Quintessence 

/ / 

/ Dust / 

m / 

Radiation 




aft) 



Phantom Matter j 
/ 



/ 



t 5 



0,2 



0.4 



0,6 0,8 



1 t 



FIG. 9: Evolution of the Universe filled with per- 
fect fluid and dark energy. 



FIG. 10: Evolution of the Universe filled with 
phantom matter. 



As one sees, equation (13.591) imposes no restriction on the value of W. But it is not the case, 
when one solves (I3.55bl) . Indeed, inserting r o ° from (12. 15al) into (I3.55bl ) one finds 



a = (A\t + Ci 



a/3(i+w) 



(3.60) 



where A\ = (1 +W) ^J3kvCq +w /4 and C\ = 3(1 + W)C/2 with C being some arbitrary constant. 
This solution identically satisfies the equation (|3.55al) . As one sees, case with W = — 1, cannot 
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be realized here. In that case one has to solve the equation (|3.55bl) straight forward. As far as 
phantom matter (W < —1) is concerned, there occurs some restriction on the value of C, as in this 
case A i is negative and for the C\ to be positive, C should be negative. As one can easily verify, in 
case of cosmological constant with W = — 1 Eqn. (I3.55bl) gives 



a e 



±^KV/3t 



(3.61) 

Inserting (|2.20al ) and (|2.20b| ) into (13.571 ) in case of Chaplygin gas we have the following equa- 



tion 



k 2Aa^+y)-XC l + y 



1+y \ 7/(1+7)- 



(3.62) 



We solve this equation numerically. The corresponding solution has been illustr ated in Fig. [£ ] 

Fin ally we consider the case with modified quintessence. Inserting (|2.24ab and (|2.24bl ) into 
(13.571) in this case we find 



K 

d =6 



(3W-1)tjCq a 



\-W JSW-2 



2W 
\-W~ 



-.e„a 



with he solution 



d 2 



vCl~ w a m - 1 



W 



l-W 



£ cr a 2 +£ 2 



E% = const. 



(3.63) 



(3.64) 



It can be shown that in case of modified quintessence the pressure is sign alternating. As a result 
we have a cyclic mode of evolution. 



IV. SINGULARITY PROBLEM 



On the the main problems of modern cosmology is the singularity problem. Let us study this 
problem within the scope of models discussed above. As we see, the components of the spinor 
field and metric functions are expressed in terms of v. It should be noted that v plays one of 
the central roles in studying the singular space-time points of BI cosmological models. Here we 
describe it in brief. In doing so let us first write the Kretschmann scalar for the Bianchi type- VI 
metric. Taking into account that the metric (13.11 ) possesses the following non-trivial components 
of Riemann tensor: 



ftm — 



R 



R 



01 

oi — 


a\ 
a\ 


p 02 

K 02 - 


a 2 


12 

12 — 


mn 
a 2 


d\ a.2 
a\ a-i 


r,13 

K 13 


31 

10 _ 


m /ci\ 
al Vai 


--)• 


r>01 

K 13 


32 

20 _ 


n /a 3 
al \ci3 


--)• 


r,02 
^23 



n03 

^03 
m 



a 



m 



«3 

3 «3 «1 

«3 a\ 



a2 a-i 



23 



a 2 ai 



(a 2 ai 
n[ 

\a2 ai 



for the Kretschmann scalar we find 
X = R„ v „rR» V(xP = 









'd\ 


4 


[( 






\a\ 


4 




(m 




a\ 





+ 



02 

a 2 

2„2 



+ 



«3 
«3 



+ 



a\ a2 
a\ a2 

2\;2 



+ 



a2 ai 
a2 ai 
d\ 

/ 

a\ 



+ 



ai a\ 
ai a\ 



a2 

«2 



(4.1) 
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The metric functions for BVI, BVIo, BV and Bill can be expressed as 

?,V /Si . (4.2) 

which gives 

(4.3a) 

h-l){y) • (43b) 

The metric functions and their derivatives for the BI take the following form: 

- = — + -^, (4.4a) 
at 3V 3V 





n ■ — 


d i _ < 






ly, 


di 




_ — t 




ai 


'V 



Clj 



, _ V 2/V\2 Xi V Xf 



a 



W~9\v) "9^ + #- ,4 ' 4W 



We study the singularity on the basis of Kretschmann scalar and in doing so we follow the 
criteria given in (Bronnikov et al, 20 04): 

(i) For any finite t some a, — > 0. (14.11) shows t hat i f more than one scale factor becomes trivial 
at finite t, then it is a singularity. As is seen from (14.41 ) as V — > (i) all a, — >■ if X\ = X2 — X3 = 0, 
i.e., it is a singularity; (ii) more than one a,- — » if more than one Xj < and in this case we have 
singularity; (iii) only one a,- — > if only one Xi < 0, i.e., the space-time can be non-singular. 

Note that this criteria of singularity always fulfills at the point where V = 0. 

As far as t — >■ °° is concerned, the corresponding asymptote can be singular if at least one a ; 
vanis hes f aster than exponentially. As X\ +X2 +X3 = 0, it means one or more Xi is negative, and 
from (14.31 ) it follows that at least one function at vanishes faster than exponentially. Hence it is a 
singularity. 

As far as BVI, BVIo, BV, Bill and FRW models are concerned, they are all singular at a space- 
time point, where V = 0. Moreover, in cases of BVI, BVIo, BV and Bill at least two scale factor 
a,'s are directly related, hence there is always a possibility of more than one one scale factor 
becoming trivial, thus giving rise to a singularity at finite t. 

Moreover, components of the spinor field as well as the physically observable quantities such as 
charge, current, spin etc. constructed from them are inverse functions of V (Saha, 2001a). Hence 
we conclude that within the scope of the model considered here at any point where V = there 
occurs a space-time singularity. 



V. PROBLEM OF ISOTROPIZATION 

Since the present-day Universe is surprisingly isotropic, it is important to see whether our 
anisotropic BI model evolves into an isotropic FRW model. Isotropization means that at large 
physical times t, when the volume factor V tends to infinity, the three scale factors ai(t) grow at 
the same rate. Two wide-spread definition of isotropization read 

E 2 = i^H 2 ^0. (5.1b) 

Here stf and E 2 are the average anisototropy and shear, respectively. Hi = di/ai is the directional 
Hubble parameter and H = d/a average Hubble parameter, where a{t) = V 1 ^ is the average scale 
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factor. Here we exploit the isotropization condition proposed proposed in Bronnikov et al (2004): 

— > const. (5.2) 



a 



Then by rescaling some of the coordinates, we can make a,- /a — > 1, and the metric will become 
manifes tly i sotropic at large t. It can be shown that in case of BVI, BVIo, BV and Bill models the 
criteria (15.21 ) does not hold, hence the isotropization process for these models does not take place. 
So we co nsider the BI model and study the problem od isotropization for this model in details. 
From (l347T) we find 

As is seen from (13.471 ) in our case a ; /a — > Z), = const as v — > °°. Recall that the isotropic FRW 
model has same scale factor in all three directions, i.e., a\(t) = ci2(t) = as(t) = a(t). So for the 
BI universe to evolve into a FRW one the constants D,'s are likely to be identical, i.e., D\ = D2 = 
£>3 = 1 . Moreover, the isotropic nature of the present Universe leads to the fact that the three other 
constants X, should be close to zero as well, i.e., << 1,0' = 1,2,3), so that Xi J [v(t)]~dt — > 
for t < 00 (for V(t) = t n with n > 1 the integral tends to zero as t — > °° for any Xj). It can be 
concluded that the spinor field Lagrangian with W < 1 leads to the isotropization of the Universe 
as t — > 00, moreover, in case of W < the system undergoes an earlier isotropization. 



VI. DISCUSSION 

Let us now examine what kind of advantage one gets exploiting the spinor description of matter. 
We do it within the scope of a BI cosmological model. In doing so, we consider the case when the 
Universe is filled with, say, Van-der-Waals fluid, radiation and quintessence. In this case we have 

r ° = e v + e r + e q , (6.1a) 

7? = -Pv-Pr-Pq- (6.1b) 

To solve (13.271) one has to know r o ° and in terms of v. It can be done exploiting Bianchi identity 

V 



T 0+y[ T l 



°-r 1 1 )=0. (6.2) 



The dark energy is supposed to interact with itself only, so it is minimally coupled to the grav- 
itational field. As a result, the evolution equation for the ene rgy de nsity decouples from that of 
the perfect fluid. Taking this into account and inserting (16.11) into (16.21) , we obtain two balance 
equations 

+ 77 fe v + Py) + £r + 77 ( £r + Pi) =0, (6.3a) 



V \ 1 J 1 V 

e q + ^(e q + Pq) = 0. (6.3b) 

In usu al approach we are in trouble, as neither e v nor e r cannot be expressed in terms of v from 
(I6.3al) . But if one uses spinor description for radiation, thanks to spinor field equation one finds 

er+-(er+p^ =0. (6.4) 
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As a result, in place of (I6.3al ) we now have 



V 

£ v + — ( £v + Pw 



0, 



(6.5) 



which is quite computable for the given equation of state. For Van-der-Waals fluid we use the 
following equation of state 

j £v 

where W\ is a constant. The dark energy density £ q can be expressed in terms of V either solving 
(I6.3bl) or using the spinor description. Thus using the spinor description (I2.15al ) and (|2.15bl ) for 
radiation and dark energy, and defining the Hubble parameter, we find the following system of 
equations 



V = 3HV, 



K 



H = -3/T + - 



8Wi£ v 



+3£ v 2 + ^y- 4 / 3 + (i-w)£ q0 y- 1 - 



£ v + 



2 

8Wl£y 

3-£ v 



w 



3ef, 



H. 



(6.7a) 
(6.7b) 

(6.7c) 



For simplicity, we set £ro = 1 and £ q o = 1 . The equation (I6.7cl) can be exactly solved for W\ = 1 /2. 
In this case we find 

£y 6 (3£v-7) Co 

As one sees, it is quite a complicate expression. In what follows, we solve the system (16.71) 
numerically. For this purpose we also set W — — 1/2 for quintessence. 
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FIG. 11: Evolution of the V, H and £ v and p for FIG. 12: Evolution of the V, H, e v and p for 

V (0) = 1, H(0) = 0.0000945 and £ v = 0.957 V(0) = 1, H(0) = -0.9 and £ v = 0.7 

In the figures [TT] and [12] we plot the evolution of V, H and £ v and p for different initial values. 
As one sees, the pressure has negative value at the initial stage of evolution, then it begin to 
increase, thus giving rise to deceleration and finally becomes negative that results in the late time 
acceleration of the Universe. 



Nonlinear Spinor Fields and its role in Cosmology 



21 



VII. CONCLUSION 

Within the framework of cosmological gravitational field equivalence between the perfect fluid 
(and dark energy) and nonlinear spinor field has been established. It is shown that different types 
of dark energy can be simulated by means of a nonlinear spinor field. Using the new description of 
perfect fluid or dark energy evolution of the Universe has been studied within the scope of a BVI, 
BVIo, BV, Bill and BI anisotropic models as well as isotropic FRW model. The corresponding 
Einstein equations have been solved. It is shown that all the models give rise to a space-time singu- 
larity where V is trivial. Among the Bianchi models considered only the BI allows isotropization 
of the initially anisotropic space-time. It is shown that the spinor description of fluid allows one to 
solve the two-fluid system without without imposing any additional condition. 
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